Rotational properties of non-dipolar and dipolar Bose-Einstein condensates confined 

in annular potentials 
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We investigate the rotational response of both non-dipolar and dipolar Bose-Einstein condensates 
confined in an annular potential. For the non-dipolar case we identify certain critical rotational 
frequencies associated with the formation of vortices. For the dipolar case, assuming that the 
dipoles are aligned along some arbitrary and tunable direction, we study the same problem as a 
function of the orientation angle of the dipole moment of the atoms. 
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I. INTRODUCTION 

In recent experiments on ultra-cold atomic gases the 
confinement of atoms in toroidal traps has become pos- 
sible [l[ . In such confining potentials persistent currents 
have also been created and observed 0t3], m close anal- 
ogy to semiconductor nanostructures 0]- The sim- 
plicity of these systems makes them ideal for studying 
effects associated with superfluidity, including rotational 
properties, persistent flow, etc. 

In another series of recent experiments it has be- 
come possible to trap atoms [§| and molecules [9| with 
a nonzero dipole moment. Contrary to most of the pre- 
viously realized experiments, in which the particles had 
zero dipole moment and the usual contact interaction 
modeled the effective atom-atom/molecule-molecule in- 
teraction rather well, in the dipolar case the interaction 
potential is long-ranged and anisotropic. As a conse- 
quence of these differences, dipolar gases present new 
properties that allow for the study of novel and inter- 
esting effects. 

In three dimensions the head-to-tail alignment of the 
dipoles causes instabilities towards collapse [Tcj, while 
low-dimensional confinement tends to make the system 
more stable Quasi- two- or quasi-one-dimensional 

confinement thus suppresses the instabilities associated 
with the dipolar interaction and is advantageous in that 
respect. 

Recently, several properties of dipolar gases have been 
examined in effectively two-dimensional systems. Among 
these, the rotational properties of dipolar degenerate 
gases have received considerable attention. In particular, 
the rotational properties of dipolar Bose-Einstein con- 
densates have been investigated in different confining ge- 
ometries, including axially-symmetric harmonic |12| . el- 
liptic harmonic [l3j . and toroidal traps fl4T - [l8j |. Also, 
recently it has become possible to create a Bose-Einstein 
condensate of dysprosium 164, with a magnetic moment 
of roughly three times larger than chromium 52, which 
increases the dipole coupling constant by roughly a factor 



often [19|| . 

Dipolar Bose-Einstein condensates confined in toroidal 
traps have been investigated both within the mean-field 
approximation 14 1, and in the few-particle limit [2fj| . 



These systems have been shown to exhibit some very in- 
teresting effects due to the interplay between the non- 
trivial topology of the ring-like confinement and the 
anisotropic nature of the dipole-dipole interaction. More 
specifically, when the dipole moments have a non-zero 
in-plane component, the effects associated with breaking 
the axial symmetry of the system have been shown to 
have potential applications in, e.g., Josephso n-ty pe oscil- 
lations and self-trapping phenomena [lg, LLZl ■ 

Motivated by the experimental progress mentioned 
above on toroidal traps and on dipolar gases, we study 
below the rotational properties of a Bose-Einstein con- 
densate that is confined in an annular-like potential, first 
considering a condensate that consists of non-dipolar 
atoms [22H24j | and then addressing the dipolar case. Our 
calculations are closely- related with those of Ref. [2o| . 
which analysed the lowest-energy state of a uniform (non- 
dipolar) rotating superfluid confined in a hard-wall an- 
nulus. In this study it was shown that as the rotational 
frequency increases initially there are vortices that are 
located in the region of zero density. As the rotational 
frequency increases, vortices start to form in the region 
of nonzero density. With increasing angular frequency of 
the trap the real vortices eventually form a circular ar- 
ray, or even multiple arrays of vortices. These theoretical 
predictions were in good agreement with the experimen- 
tal results on liquid Helium that followed afterwards (2|| . 
In what follows below we refer to the vortices that are 
located in the region of negligible density as "phantom" 
vortices and the ones that are located in the region of 
non-negligible density as "real" vortices. 

The paper is organized as follows. In Sec. II. we present 
our model. In Sec. Ill we investigate the rotational re- 
sponse of a non-dipolar gas. In the spirit of Ref. [25[ we 
identify the corresponding critical rotational frequencies. 
In Sec. IV we investigate the effect of the dipolar inter- 
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FIG. 1: (Color online) Schematic illustration of the rotating 
quasi- two-dimensional annular trap. The atoms move on the 
x-y plane with an angular frequency Q, which is along the 
z axis, in an annulus of mean radius R and a width that is 
set by the oscillator length ao = y'fi/fMwo). In the case 
of dipolar atoms, their dipole moment is oriented on the x-z 
plane, forming an angle O with the x axis due to an external 
(magnetic or electric) field E. 



action, assuming that the dipole moment of the atoms 
is oriented along some arbitrary direction due to an ex- 
ternal polarizing field. We demonstrate that the vortex 
structures depend strongly on the orientation angle of 
the dipole moment of the atoms. Finally in Sec. V we 
summarize our results. 



II. MODEL 

The annular potential that we consider is modeled via 
a combination of harmonic potentials in the transverse 
direction and in the direction perpendicular to the plane 
of motion of the atoms, 
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V(r) = V r (T ± ) + V z (z) = ~M<4(r± - R) 2 + ^M^z 



Here z is the symmetry axis of the trap, which is also the 
axis of rotation of the annulus. In addition r j_ is the posi- 
tion vector on the x-y plane, while loq and uj z are the trap 
frequencies. The ratio w 2 /wo is chosen to be equal to 100, 
which makes the motion of atoms quasi-two-dimensional, 
since Tiuj z is the largest energy scale in the problem. The 
ratio R/ao is chosen equal to 4, where ao = \Jh/ {Mlo ) 
is the oscillator length that corresponds to the frequency 
wo and mass M. Figure 1 shows schematically the corre- 
sponding annular-like potential and also the external field 
that is assumed to polarize the (dipolar) atoms. Without 
loss of generality, this field is taken to be on the xz-plane, 
forming an angle with the a;- axis. 

Concerning the interactions, these are modeled via a 
short-ranged contact potential - which takes care of the 
short-range correlations - and the usual dipole-dipole po- 
tential, which describes the long-range part of the inter- 
action. Because of the assumed strong confinement along 



the z-axis, one may safely assume that the order param- 
eter has a product form 



*(r) = $(rj_) x <j>{z), 



(2) 



where <f)(z) is the Gaussian ground state of the harmonic 
potential V z (z). This assumption allows us to reduce 
the problem from three dimensions to two dimensions 
for the order parameter $(rj_), which satisfies the Gross- 
Pitaevskii-like equation 



2M 



+ V r (rx) + Klip(rx) + glHrj-W 



*(rx) - 
= M*(r±)-(3) 



In the above equation g = (h 2 /M)y/8TrNa/a z , where N 
is the atom number, a is the scattering length for zero- 
energy elastic atom-atom collisions, a z is the oscillator 
length in the z direction, fi is the chemical potential, and 
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V eS (r±~r'±) mr' ± )\ 2 dr' ± 



(4) 



is the effective dipolar interaction potential, where 
V e s(r± — r'j_) is given in Ref. [27j . 

The problem thus reduces to solving the nonlocal and 
nonlinear integro-differential Eq. ©. We solve it with 
use of a fourth-order split-step Fourier method within 
an imaginary-time propagation approach 28[. Within 



this method one starts with some initial state, which 
then propagates in imaginary time until numerical con- 
vergence is achieved. 



III. NON-DIPOLAR ATOMS IN A ROTATING 
ANNULAR TRAP 

We first consider the case of non-dipolar atoms and 
vary the rotational frequency f2 for a fixed interaction 



(1) strength, choosing Mg/Ti — \/8irNa/a z = 150. Since 



the ratio between the interaction energy and huio is 
roughly Naao/(a z R) ~ 10, and also uj z /ljq — 100, there 
is a clear hierarchy of energy scales, with Tiui z being 
roughly ten times as large as the interaction energy, 
which in turn is roughly ten times as large as Hujq . There- 
fore, while the cloud is in the ground state of the har- 
monic oscillator along the z direction, it is closer to the 
Thomas-Fermi limit in the direction along the plane of 
motion. 

Using the numerical procedure outlined above, the en- 
ergetically favourable state of the cloud is determined 
in the rotating frame for a certain value of fl. As in 
Ref. [251 ] we identify three critical rotational frequencies 
^ci, ^c2, and C 3; see Fig. 2. The first critical frequency, 
fi c i, corresponds to the formation of a "phantom" vor- 
tex state, which is located at the center of the annulus. 
These "phantom" vortices are effectively states of quan- 
tized circulation with integer quantum number. In this 
case the density of the cloud remains unchanged and the 
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FIG. 2: (Color online) Density and phase of the order pa- 
rameter corresponding to the three critical frequencies for the 
non-dipolar case (D — 0) and \/S%Na/a z — 150. For (a) 
fl/wo = n c i/o;o ~ 0.036 a "phantom" vortex appears, for 
(b) fi/cJo = Q C 2/too « 0.53 the first "real" vortex enters the 
cloud, and for (c) Q/ljo = Q. c z/ujo ~ 0.63 an array of vortices 
forms along the annulus. 

only change takes place in the phase of the order param- 
eter, which acquires a jump of 2-7T. 

The second critical frequency, f2 C 2, is associated with 
the formation of the first "real" vortex. According to our 
results, this vortex state actually penetrates the cloud 
from the inner circle. The point where the vortex enters 
the cloud is arbitrary, as the symmetry of the Hamil- 
tonian is broken spontaneously. Eventually more vor- 
tices penetrate the cloud and finally at the third critical 
frequency f2 C 3 an array of vortices forms. We identify 
r2 C 3 as the frequency for which the vortices are located 
(roughly) at the middle of the annulus. For the parame- 
ters given above f2 c i/cjo, f^/wo, and fi C 3/wo are found 
to be ~ 0.036, 0.53, and 0.63, respectively. 

A direct comparison of the above numbers with those 
of Ref. 25] is not possible since in that case the annular 
potential was modeled as a hard-wall potential, while in 
our study we have assumed harmonic confinement and 
as a result the density varies in the transverse direction. 
According to [25j , neglecting logarithmic factors of order 
unity, r2 c i/ajo is ~ a^/(2Rd) : where d is the width of the 
annulus. Assuming that d = ao, then Q c i/ujq = 1/8, 
which is roughly a factor of 3.5 larger than the value 
found numerically. Also, according to the same study 
^2/^0 = 2(arj/d), which is equal to two, and this differs 
by a factor of roughly four from the value found above. 
On the other hand, the ratio Q, c i/Vt C 2 = ao/(4R) = 1/16, 
which is rather close to the ratio 0.036/0.53 « 1/14.7. 



IV. DIPOLAR CASE 

We turn now to the case where the atoms have a 
nonzero dipole moment. In addition to the hard-core po- 
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FIG. 3: (Color online) Density and phase of the order param- 
eter corresponding to the critical frequency fl C 2 for the dipolar 
case, with y/87vNa/a z = 150, and a^d/ao = 2.505 2 /3. For (a) 
9 = 0°, n c2 /wo « 0.80, (b) 6 = 30°, n c2 /wo « 0.5835, (c) 
9 = 60°, n c2 /cJo w 0.36, and (d) 9 = 90°, n c2 /uj « 0.296. 



tential considered in the previous section - which is kept 
fixed, with Mg/h 2 = V^Na/a z = 150 - we consider a 
dipolar interaction of a fixed strength. Introducing the 
dipolar length add = MD 2 /(3K ), we choose add/a = 
2.505 2 /3. In the above expression D 2 = d 2 /(4ireo) when 
the atoms have an electric dipole moment d, where eo 
is the permittivity of the vacuum; when the atoms have 
a magnetic moment /i, D 2 = /io/i 2 /(47r), where fig is 
the permeability of the vacuum. With this choice of pa- 
rameters the dipolar and the contact interactions are of 
comparable magnitude. 

In this problem, small values of D and thus weak dipo- 
lar interactions may give rise to a very smooth energy 
surface, which leads to degeneracy problems in the cal- 
culation. On the other hand, large values of D tend to 
cure this problem, however they also make the system 
unstable against collapse for small values of the angle O 
due to the head-to-tail alignment of the dipoles. 
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We investigate the rotational response of the cloud con- 
sidering four values for the angle of the dipole moment, 
9 = 0°,30°,60°, and 90°, as in Fig.l. When 6 = 90° 
the interaction respects the axial symmetry of the trap- 
ping potential and is also purely repulsive. As a result, 
the system behaves qualitatively as in the case of a con- 
tact potential alone, with a phantom vortex state, a real 
vortex state, and a vortex lattice forming with increasing 

n. 

For any other value of the angle 9, as fl increases, 
initially the cloud still responds by forming one phantom 
vortex at the center of the trap, as in the non-dipolar 
case. On the other hand, for higher rotational frequencies 
the gas behaves in a qualitatively different way. In this 
case the dipolar interaction breaks the axial symmetry of 
the Hamiltonian, introducing a preferred direction. 

In Fig. 3 we plot the density and the phase of the order 
parameter for four values of the angle 9 = 0°,30°,60°, 
and 90°, for the value of the angular frequency that cor- 
responds to f2 C 2, i.e., for the value of ft for which the 
first (real) vortex state penetrates the annulus. While 
for 9 = 90° the position where the vortex enters the 
annulus is arbitrary (as in the non-dipolar case), for all 
other values of O this is determined by the direction of 
the polarizing field, which is chosen to be in the direc- 
tion going from the bottom to the top of the page. For 
9 = 60° and 30° there are actually two vortex states 
which enter the annulus from opposite ends. For 9 = 0° 
instead of two vortices, there are two pairs of vortices, 
which are symmetrically displaced from the direction of 
the polarizing field. As one can see from the phase of the 
order parameter, as decreases the number of phantom 
vortices increases and especially for 9 = 0° their number 
is very large. Also, with decreasing 9 the density at the 
center of the trap decreases. For 9 = 0° the density in 
this region is extremely small and therefore it appears 
to be "empty" in Fig. 3, as it can be seen also from the 
phase plot. 

Furthermore, as we have seen numerically, the 
anisotropic character of the dipolar interaction affects 
the frequencies O c i and Q C 2 in a different way. For 
9 = 0°,30°,60° and 90°, Q el is found to increase lin- 
early with 9, with Q c i/ujq being approximately equal to 
0.035, 0.036, 0.037, and 0.038, respectively. Also, fl c2 /u} 
is found to be « 0.80, 0.5835, 0.36, and 0.296, respec- 
tively. In other words, as 9 decreases, i.e., as the system 
gets more distorted from axial symmetry, Q c i decreases 
slightly (and linearly), whereas fi C 2 increases. It is inter- 
esting that for 9 = 30°, fl c i is identical to that of the 
non-dipolar gas. For this value of 9 the repulsive and at- 
tractive parts of the dipolar interaction are comparable 
and in a sense balance each other. 



V. SUMMARY 

In the present study we have investigated two prob- 
lems. First, we studied the rotational response of a su- 
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FIG. 4: (Color online) The critical frequencies il C 2 and £l C 3 
as a function of the angle 6 for the parameters of Fig. 3. 

perfluid that is confined in an annular potential. Such a 
trapping geometry interpolates in a sense between one- 
and two-dimensional motion, depending on the ratio be- 
tween the mean radius of the annulus and its width. The 
second problem we have addressed is the effect of the 
long-range and anisotropic character of the dipolar inter- 
action. 

Starting with the non-dipolar the rotational 

frequency of the annulus increases, there is a critical fre- 
quency above which a vortex state appears in the region 
of essentially zero density, around the center of the trap. 
Eventually more vortices start to reside in this region. 
These vortices affect only the phase of the order param- 
eter, increasing the circulation and leaving the density 
unaffected. As the rotational frequency increases, the 
number of "phantom" vortices increases as well. Even- 
tually the number of phantom vortices which can be ac- 
commodated in the region of low atom density reaches a 
threshold. As a result, above some critical frequency, one 
of these vortices penetrates the annulus, moving to the 
region of non-zero density. These "real" vortices eventu- 
ally form a regular vortex array, or even multiple arrays, 
with increasing ft. 

The dipolar case differs from the above mainly because 
of the symmetry-breaking of the Hamiltonian due to the 
dipolar interaction (unless the dipole moment is perpen- 
dicular to the plane of motion of the atoms). Even the 
non-rotating system is affected by the dipole-dipole in- 
teraction, with the formation of one or two density max- 
ima, depending on the value of the parameters, as seen 
in Refs. [H, [3 EH- This inhomogeneity in the density 
affects also the critical rotational frequencies. Initially a 
"phantom" vortex still appears at the center of the trap. 
However, for larger rotational frequencies of the trap the 
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"real" vortices penetrate the annulus either in pairs, or 
in couples of pairs from opposite ends of the cloud in the 
direction of the polarizing field of the dipoles. 

Rotating atomic Bose-Einstein condensates confined in 
toroidal potentials have been realized experimentally re- 
cently, while dipolar gases are under intense experimental 
investigation. As shown in the present article, the com- 
bination of these two effects - that should be accessible 
also experimentally - gives rise to interesting effects. 
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